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ABSTRACT 



Aims. We model a diffuse molecular cloud present along the line of sight to the star HD 102065. We compare our 
modeling with observations to test our understanding of physical conditions and chemistry in diffuse molecular clouds. 
Methods. We analyze an extensive set of spectroscopic observations which characterize the diffuse molecular cloud 
observed toward HD 102065. Absorption observations provide the extinction curve, H2, C I , CO, CH, and CH + column 
densities and excitation. These data are complemented by observations of C + , CO and dust emission. Physical conditions 
are determined using the Meudon PDR model of UV illuminated gas. 

Results. We find that all observational results, except column densities of CH, CH + and H2 in its excited (J > 2) levels, 
are consistent with a cloud model implying a Galactic radiation field (G ~ 0.4 in Draine's unit), a density of 80cm -3 
and a temperature (60-80 K) set by the equilibrium between heating and cooling processes. To account for excited 
(J > 2) H2 levels column densities, an additional component of warm (~ 250 K) and dense (hh > 10 4 cm -3 ) gas within 
0.03 pc of the star would be required. This solution reproduces the observations only if the ortho-to-para H2 ratio at 
formation is ~ 1. In view of the extreme physical conditions and the unsupported requirement on the ortho-to-para 
ratio, we conclude that H2 excitation is most likely to be accounted for by the presence of warm molecular gas within 
the diffuse cloud heated by the local dissipation of turbulent kinetic energy. This warm H2 is required to account for 
the CH + column density. It could also contribute to the CH abundance and explain the inhomogeneity of the CO 
abundance indicated by the comparison of absorption and emission spectra. 
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1. Introduction 



shock s dDraine fc Katl Il986t 



Since the pioneering work of iBlack fc Dalgarnol (|l977t ). ob- 
servations of diffuse molecular clouds continue to motivate 
and challenge efforts to model the thermal balance and 
chemistry of interstellar gas illuminated by UV photons. 
Models allow observers to determine physical conditions 
from their data and observations contribute to models by 
quantifying physical processes of general relevance to stud- 
ies of matter in space such as H2 formation, photo-electric 
heating, and cosmic ray ionization. 

Many models of well characteriz ed lines of sight have 
been presented (e.g. in the last years: Zsargo fc Federman I . 
120031 : iLe Petit et all . 120041 : IShaw et all . I2006T ). Thev are 
successful in reproducing many observables apart from 
some molecular abundances, most conspicuously CH + , 
which points to out-of-equilibrium chemistry. This molec- 
ular ion, and several of the molecular species com- 
monly observed in diffuse molecular clouds such as 
CH, OH and HCO+ may be produced by MHD 



Pineau des Forets et all 



119861 : iFlower fc Pineau des Forets . 1998H. and sma ll scale 
vortices ~(| Joulain et all . 1998^ Falgarone et ahl . l2006f) where 
H 2 is heated by the localized dissipation of the gas turbu- 
lent kinetic energy. Turbulent transport between the cold 
and warm neutral medium ma y also significantly i mpact 
the chemistry of diffuse clouds (jLesaffre et al.1 . 120071 ) . 

Independently of gas chemistry, the presence of H2 
at higher temperatures than that set by UV and 
cosmic-rays heating of diffuse molecular clouds, may be 
probed through observations of th e H 2 level popula- 
tions (jCecchi-Pestellini et al.1 . 120061 ). A correlation be- 
tween CH + and r otatiq nally excited H 2 was found by 
Lambert fc Dankd (1 1986ft using Copernicus observations. 



Falgarone et al.l (|2005ft reported the detection of the S(0) 
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to S(3) H 2 lines in a line of sight towards the inner 
Galaxy away from star forming regions. They interpret 
their observation as evidence for traces of warm molecu- 
lar gas in the diffuse interstellar medium. But the inter- 
pretation of the wealth of H 2 observations pr ovided by 
the F USE satellite is still a matter of debate. I Gry et all 
(|2002ft modeled FUSE H 2 observations of three stars in 



Chamaeleon usi ng the Meudon Photon D ominated Regions 
(PDR) model (|Le Bourlot et all Il993f ). They show that 
the model cannot account for H2 column densities in ro- 
tational stat es with J > 2. A l arger sample of H2 FUSE 
observat i ons (iTumlinson et al. I , l2002t iGillmon et "all , l2005t 
IWakkerl . 120061). inclu d ing 2 of the 3 Chamaeleon lines 
of sight of iGrv et al.1 (|2002t ). have been analyzed on the 
basis of model calculations presented by iBrowning et al.l 
(|2003f ). Their model, like other PDR models, takes into ac- 
count the formation of H2 on grains, its photo-dissociation 
by absorption of resonant UV photons, radiative transfer 
and vibrational/rotational excitation resulting from colli- 
sions, H2 for mation and UV pumping. Unlike the l atest 
PDR models (IShaw et all 120051 ; iLe Petit et ail . l2006f ) , the 
IBrowning et all (|2003l ) model docs not derive the gas tem- 
perature from th e thermal balance betw een heating and 
cooling processes. IBrowning et al.l (|2003l ) instead consider 
gas temperature as a model parameter independent of the 
density, incident radiation field and cloud shielding (total 
extinction). They conclude that H2 observations cannot be 
accounted for with a single isothermal slab of gas. They 
propose solutions where the data are fitted with absorption 
from two physically independent gas layers with distinct 
temperatures, UV illuminations and column densities. A 
cold (< 100 K) component contributes most of the total 
H2 column density while a warmer (~ 200 K) and thin- 
ner component with a higher UV field helps populate the 
high J s tates. For many of the model combinations consid- 
ered by IBrowning et all (|2003l ). the difference in UV field 
is too low to account for the corresponding difference in 
temperature. In particular, the combinations proposed for 
the nearby diffuse interstellar medium require an additional 
heating source in the warm component. For H2 observa- 
tions towards Galactic stars, some of the ab sorption may 
arise from gas in the vicinity of the star (e.g. iBoisse et all 
120051) . 



In a companion paper (jNehme et all l2008h . we pre- 
sented a multi-wavelength study combining spectroscopic 
UV, optical, IR and radio observations of the interstel- 
lar matter, along the line of sight to the nearby (170 pc) 
moderately reddened (E(B-V) = 0.17) star HD 102065. 
Absorption observations provide the column densities of H2 
in the J=0 to 5 states, C 1 in its three fine structure states, 
CO in the J=0 to 2 states, CH and CH + . They are comple- 
mented by observations of the C + , CO(2 — 1) and (1 — 0) 
lines and dust continuum in emission. Non-detections of the 
C2 and CN optical absorption lines were transformed into 
upper li mits on column densities u sing oscillator strengths 
listed bv lGredel etaLl (j!991l Il993h . 

The HD 102065 line of sight is well suited for de- 
tailed modeling of physical conditions, chemistry and H2 
excitation, because of the large amount of available data 
(Table [T]). The molecular fraction is not altered by the pres- 
ence of warm atomic gas along the line of sight. Comparison 
of H I, H2 and dust extinction indicates that the bulk ( 90%) 
of the column density is accounted for by a diffuse molec- 
ular cloud identified on IRAS images. IRAS observations 
place a strong constraint on the presence of matter close 
to the star. The a bundance of CH + is high. This paper ex- 
tends the work of IGrv et all (|2002t ). where only the FUSE 
H2 spectrum of HD 102065 was analyzed, to a wider set of 
observ ations, using an updat e d version of the Meu don PDR 
model (|Le Petit et al.l . l2006h . iKopp et al. I (|2000h used this 



sight line to discuss the impact of far-UV extinction on the 
CO abundance. 

The structure of the paper is as follow: Sect. [3] presents 
the PDR model used to characterize the diffuse molecu- 
lar cloud and Sect. [3] describes the main modeling results. 
First, a reference model is defined by fitting observational 
constraints together. Second, the model is compared with 
each of the observations to assess the dependence of model 
predictions on the values of the physical parameters. In 
Sect. HI we present a detailed attempt to model H2 exci- 
tation with a warmer component located close to the star. 
Sect. [5] presents our conclusions. 



2. Description of the PDR model 

We use a comprehensive model of an interstellar cloud, that 
describes the state of the gas and dust exposed to a ra- 
diation field as a function of optical depth . The model 
is one-dimensional and stationary. It pr esents seve r al im- 
provements over that previo usly used by IGrv et al.l (|2002t ) 
and is described in detail in ILe Petit et al.l (|2006l ). It com- 
putes simultaneously, in an iterative way : 

— the UV radiative transfer, taking into account the 
continuum absorption by dust and the line absorption 
by the lowest levels of H2. The UV radiation field can 
illuminate the cloud from both sides. 

— the thermal balance, taking into account all relevant 
heating and cooling processes. 

— the chemistry, typically coupling about 100 different 
species through 800 chemical reactions. 



2.1. Model parameters and variables 

Model parameters (see Table[5|) are kept fixed at values con- 
sistent with typical diffuse clouds and the measured char- 
acteristics of the HD 102065 line of sight. The model ig- 
nores the gas velocity structure and as sumes that th e three 
low velocity components discussed by iNehme et al.l (|2008[) 
make a single homogeneous cloud with a visual extinction 
equal to the observed value to HD 102065. We look for 
a best fit model by varying only the UV radiation field 
strength G (isotropic and incident on both sides of the 
cloud) measured in units of Draine's radiation field, and the 
gas density nn- Density is constant throughout the cloud. 
Temperature is computed by solvi ng for therma l balance. 
The standard Draine's field (G = l. lDraind . )l978l) is equiv- 



alent to about 1.6 in units of Habing's field ( Habind . ll968l ). 

Element abundances are t h ose m easured for £ Oph 
taken from ISavage fc Sembachl (|l996h . Fig. [1] shows that 
HD 102065 measured abundances are close to these refer- 
ence values (Paper I). The grain distribution is kept fixed. 
It is a MRN type one with a power law size distribu- 
tion a = —3.5 with minimum and maximum sizes 10~ 7 
and 1.5 x 10 _5 cm. Although that distribution does not 
reproduce the details of the smallest grains size distribu- 
tion required to account for the mid- i nfrare d dust emission 
(|Desert et all Il990t iBoulanger et all . Il994l ) , it matches the 
Very Small Grains (VSGQ) to total dust mass ratio of 0.25. 



1 By VSG we mean any solid particle with a typical size 
smaller than 10 -6 cm, irrelevant of its precise nature. 



Table 1. Observational constraints and best model results. 
Upper part are constraints used in Fig. [2j lower part com- 
pares unconstrained observations and results. Number in 
parentheses are powers of 10. 





j^-mod 




Cobs 


jV(CO)/jV(H 2 ) 

TV T f \ /AT 

N(C i)/Nh 
N(C i} =1 )/iV(C i) 
N(C iT =2 )/N(C i) 

./H 2 N(U) + 2N(U > ) 

AT/T TO \ /T\r/TTP\ 

I(C+) (erg/s cm 2 sr) 


1.5 (-7) 
5.8 (-7) 
0.17 
0.03 
0.9 
0.73 
2.0 (-6) 


1.6 (-7) 
6.0 (-7) 
0.2 
0.02 
0.69 
0.7 
2.8 (-6) 


,0.2 (-7) 
^O.IS (-7) 

±1.5 (-7) 
±0.07 
±0.01 
±0.12 
±0.12 

±0.85 (-6) 


iV(CN)/JV(H 2 ) 
AT(C 2 )/iV(H 2 ) 


8 4 f— 91 
1.2 (-10) 
3.6 (-8) 


1.85 (—8) 

< 1.5 (-9) 

< 3.5 (-8) 


±0.3 (—8) 


jV(COj =0 )/iV(H 2 ) 


9.0 (-8) 


9.6 (-8) 


,1.4 (-8) 
^1.7 (-8) 


AT(CO J= i)/iV(H 2 ) 
iV(CO J=2 )/iV(H 2 ) 


5.1 (-8) 
3.7 (-9) 


6.2 (-8) 
< 7.3 (-9) 


,1.5 (-8) 
^1.2 (-8) 



Table 2. Fixed model parameters. Sx are the gas phase 
abundance of atom X relative to hydrogen. 



Parameter 


Value 


Comment 




0.67 mag 


Extinction 


C 


SlO^s" 1 


Cosmic Rays ionization 


Vturb 


2.0 km s -1 


Turbulent velocity 


11 y 


3.9 


Total to selective extinction 


U! 


4.2 10 _1 


Dust albedo 


fl 


6.0 10" 1 


Dust anisotropy factor 




0.57 10" 2 


Dust to gas mass ratio 


Pg 


2.59 gem" 3 


Dust density 


O 


3.5 


Dust size distribution index 


^min 


1 10 _7 cm 


Dust minimum radius 




1.510~ 5 cm 


Dust maximum radius 


Sc 


1.32 10" 4 




So 


3.19 10" 4 




<5n 


7.50 10~ 5 




<5s 


1.86 10" 5 





X 

X 



-1 



-2 



-3 




-• — ^oph Cool 
* HD1 02065 CompA+B 



-4 

O P Zn Si Cu Mn Mg Fe Cr Ni Co Ti 
Element 

Fig. 1. Depletion [X/H] in the cool gas towards HD 102065 
compared with £ Oph values. Data are taken from Paper I. 



2.2. H2 formation equations 

In the present model, the H2 formation rate on grains is 
not a free parameter, but is computed locally from the com- 
bined adsorption of H atoms from the gas on grains followed 
by recombination and desorption of H2 molecules: 



H + dust - 

Had + H a d 



Had 

^H 2 



where H a d is an hydrogen atom adsorbed on dust. If only 
those two processes are included and steady state applies, 
the production of H2 is independent of how H atoms even- 
tually manage to reach one another and may be computed 
by: " 



rf[H 2 



dt 



1 



k ad [H] 



1 



s < crn g > vh n(H) (1) 



form 



where s is the sticking coefficient of H upon collision, Tin its 
mean velocity and < an g > the mean grain cross section 
per unit volume. 

In the full model, other processes may compete with 
H2 formation to remove adsorbed hydrogen atoms (photo 
desorption, reactions on grains, etc.). In this paper, we 
neglect all of those processes, which allo ws the formation 
rate to be computed exactly. Results from lLe Bourlot et al.l 
(| 1995T ) show that the mean cross section can be written: 



< an 9 >= 4 



3 1.4 m H G r 



Pg 



\fa~n 



x a„ 



(2) 



where G r is the dust-to-gas mass ratio, p g the density 
of grain material, and integration is completed for the 
entire grain size di stribution, using a MRN distribution 
(jMathis et al.l . [l977[ ) with an index of 3.5 from a mm to a max . 
Combining Eq (p]) and j2]) we can write: 



d[H 2 



dt 



= R f v/Tgas s n H n(H) 



(3) 



form 



where i?f, the H2 formation rate in cm 3 s , reads: 

1 



Rt = 1.27 x 10 - 20 — 



Pg V dmin x a n 



The resulting value of Ri is 2.3 x 10 



-17 



(4) 



We 



use a prescription f rom David Flower (see discussion in 
iLe Petit et al.l . [20061 ) for the sticking coefficient s: 



= ,/l0/T ga8 (K) 



(5) 



The ortho-to-para ratio on formation of H 2 is set to a 
value of three. One third (1.5 eV) of the formation energy 
is transferred into internal excitation, and is distributed 
over all H 2 energy l evels with a Boltzmann distribution 
(|Le Petit et al.l . l2006h . A second third of the formation en- 
ergy is converted into kinetic energy of H 2 , and the last 
third, into grain heating. 



3.35 1(T 



200 




Fig. 2. x 2 contours (Eq \§§ using the top 7 quantities of 
Table [TJ The best fit is obtained for G — 0.4 and njj = 
80cm~ 3 . Contours are labeled with the % 2 value. 



3. Modeling the molecular cloud toward HD102065 

3.1. General modeling results 

To compare observations and model results, we have 
selected seven quantities for which an observa- 
tional error bar could be computed: iV(CO)/iV(H2), 
N(C i)/N H , N(C ij = i)MC i), JV(C C 2 )/iV(C i), 
/h 2 = 2N(R 2 )/(N(R) + 2N(U 2 )), N(H%)/NQ%), and 
I(C + ), where N(X) is the column density of species X over 
the whole cloud and and are ortho- and para-H2 
respectively. From these quantities, wc define a x 2 error 
function by: 



x 2 = 



Cobs 



(6) 



where X° hs is the quantity derived from observations, JQ nod 
the same quantity from a model, and cr b s the observational 
uncertainty (see Table QJ 0. We computed a grid of mod- 
els for gas densities ranging from 30 to 200 cmT z and UV 
fields G from 0.25 to 1.0. Figure [2] shows the resulting x 2 
iso-values. Contours are smooth, and limit a well defined 
minimum where x 2 < 1, i.e. where model results are, in the 
mean, closer to observations than observational uncertain- 
ties. 

One can see that the best compromise is reached for a 
rather low radiation field (G = 0.4, which is close to the 
standard Habing field) and a mean total proton density 
of Tin = 80 cm -3 . In the following wc will refer to those 
values as our reference model. The reference model results 
are compared to the observables in Table [U Abundance 
and temperature profiles are illustrated on Fig.l3"ll4landlTlFl. 
The H I to H2 transition occurs close to the cloud edge at 
an extinction of A v ~ 10 -3 . Most of the carbon is in the 
form of C + at all depths. The fraction of atomic carbon 



Note that Xi need not be dimensionless because normalisa- 
tion by <r ensure that compatible quantities are summed 

3 The peculiar horizontal scale allows for logarithmic scaling 
towards both sides of the cloud 




-1 -1 

Log 10 (min(A v )/A v (c)) 

Fig. 3. H and H2 density profiles for the reference model 
plotted versus the extinction from the nearest edge normal- 
ized to the central extinction. 
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Fig. 4. C + , C, CO and CH density profiles for the reference 
model. 



C 1 is constant throughout the cloud. The density profiles 
of CO and CH molecules follow that of H2. Temperature 
varies from 60 to 80 K, and is lowest at the cloud edge. The 
temperature bumps visible on Fig. [5] are due to heating by 
H2 formation. 

The gas heating rate is 1.8 x 10 -24 erg cm" 3 s^ 1 at the 
cloud center and comes mainly from photo-electric effect 
on grains. In the outer layers where the temperature peaks, 
photo-electric heating reaches 3.2 x 10 _24 ergcm~ 3 s^ 1 and 
there is a lower but significant heating due to H2 formation 
(2.2 x 10~ 24 erg cm" 3 s" 1 ). Cooling is dominated by C + , 
with H 2 contributing up to 20% around 10 -2 A v . This leads 
to an integrated C + emissivity of 2 x 10 -6 erg cm -2 s" 1 sr~ 4 
in good agreement with the ISO observation. 

In the next sections, we discuss each of the observables 
and their dependence on the two model parameters n# and 
G. 
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Fig. 5. Temperature profile for the reference model. 
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Fig. 6. H2 fraction /h 2 as a function of njj/G for G = 
0.25, 0.4 & 1.0. The horizontal lines are the observed value 
with the error bars. 
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Fig. 7. Ortho-to-para ratio of H2 as a function of nn for 
three- values of G. 



Fig. 8. Mean gas temperature as a function of nn for three 
values of G. 



3.2. H2 formation rate and molecular fraction 

H2 fraction (fH 2 ) is plotted versus density for three val- 
ues of G on Fig. [5] For the reference model (G = 0.4, 
riH = 80 cm" 3 ) computed values for fjj 2 are around 0.9, 
significantly higher than the observed fjy 2 . 

iGrv et afl (|2002l ) have shown that /h 2 , for a given value 
of G, depends on the product nn x R, where H2 formation 
proceeds at a rate nn n(H) R. We note that their R incor- 
porates both the temperature dependence and the sticking 
coefficient. Using our prescription for the sticking coeffi- 
cient s (section ([5])), one finds that R is equal to \/T0 R{ 
and is independent of T gas . 

More generally, Fig. [6] shows that the molecular fraction 
for a given total gas column density depends on njj X R/G. 
IGrv et aT] (|2002l ) found that the observed fn 2 is matched for 
7iH X R = 2.3 10~ 15 s _1 , for a cloud illuminated on one side 
by the Habing field (equivalent to G = 0.6 in Draine units). 
Our reference model gives nn VT0 i?f = 5.8 10 _15 s _1 , 
a factor 2.5 higher. This difference explains why, in our 
model, the molecular fraction is higher than the observed 
value. In the model, the formation of molecular hydrogen is 
the only destruction path of adsorbed atomic hydrogen in 
the present calculation. Other processes (such as photo des- 
orption) may limit the amount of H adsorbed on grains and 
lower further the H2 formation efficiency. We consider the 
agreement between model and observed H2 fraction sat- 
isfactory because the time scale of H2 formation, a few 
10 7 yr, is long compared to dynamical time scales in dif- 
fuse clouds. The observed value is thus not expected to 
accurately match the model steady state value. Analysis 
of a large set of FUSE lines of sight will be required to 
estimate how observed molecular fractions scatter about 
model steady state values and constrain the H2 formation 
rate more precisely. 



3.3. H2 excitation 

The ortho-to-para H2 ratio, and the mean gas temperature 
are shown in Fig. [7] and [8j and the H2 excitation diagram 
in Fig. [Ul H2 column densities in the J = and J = 1 
levels are reasonably well reproduced. The observed ortho 
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Fig. 9. Observed H2 excitation diagram and reference 
model. Two excitation temperatures are plotted for J = 
and 1: T ex = 66 K, and for J = 3 to 5: T ox =248K. 



to para ratio depends on the gas temperature. If it is at its 
equilibrium value (this is a reasonable hypothesis for diffuse 
clouds where the time scale for ortho-to-para conversion is 
shorter than that of H 2 formation), the N(J = 1)/N(J = 0) 
ratio can be considered a direct measure of the gas temper- 
ature. The excitation diagram in Fig. [5] implies an excita- 
tion temperature of 66 K (for J< 2), which indeed is in good 
agreement with the mean gas temperature for the reference 
model (see Fig. [5]). 

On the other hand, H2 column densities derived from 
the model for J > 2 are far below the observed values. 
For levels J > 3, we derive an excitation temperature 
T cx = 248 K. This high J gas suggests the existence of a 
warmer gas component that will be discussed in Sect. [4] 
We used the model to check that an increase in the cosmic- 
rays io nization rate, a possibility considered by I Shaw et alj 
(|2006h . is not a solution. Multiplying this rate by a fac- 
tor 20 to a value of 10" 15 s _1 increases the model column 
densities in the J=3 and 4 levels by less than a factor of 2. 

3.4. Carbon abundance and excitation 

The neutral atomic carbon abundance and fine structure 
excitation, are presented in Fig. [10] and [TlJ Most of the 
carbon is in C + and the fraction of carbon in the neutral 
form C 1 depends on both the gas density and the value of 
G. It increases with densities and decreases with increas- 
ing G. Our reference model is in good agreement with the 
observed value. 

The neutral atomic carbon excitation in diffuse 
clouds is often conside red as a measure of gas pressure 
(| Jenkins fc Trippl . l200ll) . In our model, this translates into 
a dependence on the product nn x G with only small dif- 
ferences from one G to another. Our reference model, re- 
produces well the fraction of J = 1 fine structure level, but 
not so well that of J = 2, which implies n# < 50 cm -3 for 
G = 0.4. It is often the case that C 1 excitation in diffuse 
clouds cannot be accou nted for by gas at a single pressure 
(| Jenkins fc Trippl . l200lh . However our observations show a 
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n H (cm 3 ) 

Fig. 10. Neutral atomic carbon abundance as a function of 
the gas density Tin for three values of the UV radiation field 
intensity factor G. 
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CI*7Cltot 

Fig. 11. Fraction of excited C I (J=l and 2 fine structure 
levels: C I* and C I**) for various nn and G. Dashed lines 
are the observations with their 3-er uncertainties. Along 
each line, the densities increase from left to right from 30 
to 200 cm" 3 as in Fig. fTUl 



lo wer C i**/C I* rat i os tha n implied by the model, contrary 
to I Jenkins fc Trippl (|200lh findings. 

3.5. Molecules 

Our reference model matches the observed abundance of 
CO (Fig. [12]) . Tabled] shows that CO excitation is also well 
predicted by the model, although it was not used as a con- 
straint. The rise in the iV(CO)/iV(H2) ratio at low density 
and low radiation field is a chemical effect. The formation 
of CO is illustrated in Fig. [T3] The efficiency of CO for- 
mation depends on the ionization fraction. Photoionisation 
is proportional to the gas density ne, while recombination 
proceeds at a rate that is proportional to the square of that 
quantity. The degree of ionization therefore increases when 
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Fig. 12. N(CO)/N(H 2 ) variation with density. All other 
model parameters are as in Table [I] The horizontal lines 
show the observed values with the error bars . 
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Fig. 14. N(CH)/N(H2) variation with density and the ra- 
diation field. The horizontal lines show the observed values 
with the error bars. 
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Fig. 13. A schematic illustration of CO formation routes 
in a low UV radiation field. Destruction is still dominated 
by Photodissociation. 



the density decreases in a given radiation field. Enhancing 
the ionization enhances + formation via charge exchange 
with H + and favors the formation of OH and H2O. Both 
molecules interact with C + and lead either directly or indi- 
rectly to the formation of CO. The reaction involving OH is 
dominant. Photodestruction of CO is limited by the weak- 
ness of the radiation field. 

We are cautious about the model interpretation of the 
CO abundance. It is noticeable that the CO column den- 
sities in the J=l and 2 levels, derived from UV spectra, 
are significantly smaller by factors 3 and 7, respectively, 
than the valu e s der ived from the emission radio spectra 
(jNehme et al.l . l2008f ). These differences indicate that the 
abundance of CO is inhomogeneous, which could be ac- 
counted for within the PDR model, by introducing clumps 
with higher density than the mean value. However, the large 
CH + abundance favors an alternative explanation, where a 
small-scale increase in the CO abundance, traces the local- 
ized con tribution of out-of-equi librium chemistry to its for- 
mation. iLiszt and Lucasl (|2000[ ) have gathered results from 
UV and radio absorption measurements of CO along diffuse 
interstellar medium lines of sight. They relate the CO abun- 
dance to that of HCO + measurements, concluding that CO 
formation through dissociative recombination of HCO + suf- 



fice to acc ount for the CO abun dance in diffuse molecu- 
lar clouds. iFalgarone et alj (|2006T ) show that the observed 
abundance of HCO + (~2x 10~ 9 ) cannot be accounted by 
standard PDR chemistry and must be related, like CH + , 
to warm out-of-equilibrium chemistry. A significant frac- 
tion of CO observed in diffuse molecular clouds may thus 
be a product of out-of-equilibrium chemistry. 

This interpretation links CO abundance inhomo- 
geneities to the CH + chemistry but it is not specifically 
the CH + rich gas that has an enhanced CO abundance. 
It is not ruled out b y the observed velo city difference be- 
tween CH+ and CO jNehme et alll2008f ). CO, unlike CH+, 
is observed to be concentrated in the intermediate velocity 
component C. This component may correspond to shielded 
sections of th e cloud where the C O photo-dissociation rate 
is reduced. In lNehme et al.l (|2008l ). we propose that the line 
of sight to HD 102065 samples material ablated from the 
Dcld 300.2-16.9 cloud by a cloud-supernova shock interac- 
tion. In this scenario, the matter flowing out of the clo ud 
is expected to be very turbulent (jNakamura et all l2006f ) . 

In Fig. HU we show that the CH abundance 
(N(CH)/N(H 2 )) computed by the model depends linearly 
on the ratio nn/G. The value for our reference model is 
a factor of two lower than observations. This mismatch 
between model and observations, might be an additional 
manifestation of the ou t-of-equilibrium chem i stry, a s al- 
ready propos e d by e.g. IZsargo fc Federman I |2003!) and 
iRitchev et ail (|2006l ). The "excess" of CH abundance may 
be the product of CH + recombination with H 2 . Finally, we 
note that the model is consistent with the upper limits on 
CN and C 2 abundances (Table [1]). 



4. Warm H 2 

The H 2 excitation diagram (Fig.HJ shows that the observed 
column densities at J >2 lie far above the single temper- 
ature fit to the low J column densities. This is commonly 
observed, towards many stars by Copernicus and FUSE. 
The model values, which take into account UV pumping 
and H 2 excitation upon formation are also roughly an order 
of magnitude lower than the observations. iGrv et al.l (|2002h 
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Fig. 15. H2 (J = 3) column densities for an A v = 1CP 2 
slab of gas close to the star (labels next to the curves is d in 
pc). The two dashed horizontal lines bracket the measured 
column density with its error bar. 



and others earlier (e.g.lDraine fe Katzl . Il986l : iJoulain et all 



119981 : iPineau des Forets et all Il986h have proposed that 



this excited H2 traces warm gas in regions where kinetic 
energy is dissipated through shocks or vortices. Other au- 
thors have proposed that this warm H2 is the signature of 
molecular gas close to the star and thus exposed to a high 
UV field (e.g. lBrowning et alll200l . The presence of CH+ 
in quantities much larger than predicted by the PDR m odel 
favors the former explanation (jFalgarone et all l2005h . To 
test this preference, we have used the PDR model to quan- 
tify the latter possibility in the specific case of HD 102065. 

The ratio between the 100/xm brightness in the IRAS 
images and the visible extinction towards HD 102065 in- 
dicat es that most of the m atter is fairly distant from the 
star (jBoulanger et all Il994l ) and does not interact with it. 
The 60 and 100 /jm images only show a small brightness 
enhancement at the position of HD 102065, point-like at 
the IRAS resolution (angular size ^ 5'), and correspond- 
ing to a small fraction ( 10 -3 ) of the stellar luminosity 
(jBoulanger et all ll994). If the absorbing matter occupies 
a solid angle 8 about the star, its UV/visible opacity is 
10~ 3 x 4n/6. Combining the constraint on the source di- 
ameter (< 5'), and the star distance (170 pc), we derive an 
upper limit for the distance, from the star to the absorbing 
matter of 0.12 pc. HD102065 is a B9IV star with a luminos- 
ity ~ 100 Lq and an effective temperature of 11300 K. The 
stellar radiation field intensity is G* = 0.2 x (<i/lpc)~ 2 in 
Drainc units. 

For an extinction A v ~ 10~ 3 the gas is molecular only 
for large densities. To search for a more realistic solution, we 
assume that O/An ~ 0.1. To quantify with this assumption, 
the column density of warm H2 that could be associated 
with the IRAS far-IR emission, we compute a grid of models 
with constant A v = 10~ 2 , varying the density «h from 10 3 
to 3 x 10 6 cm -3 , and the distance to the star from O.lpc to 
1 10~ 2 pc (equivalent to G* from 20 to 2 x 10 3 ). 

Figure HU shows the column densities of H 2 (J = 3), as 
a function of «h for various distances to the star. As dis- 
tance decreases, the radiation field increases, the H2 photo- 
dissociation rate rises, and higher densities are required to 
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Fig. 16. H2 excitation diagram for nn = 3 10 4 cm~ 3 , d = 
0.02 pc (G ~ 300). The J=4 column density is reproduced 
only if the H2 ortho to para ratio at formation would be 1 
but no theoretical nor experimental study support a value 
much different from 3. 



retain hydrogen in its molecular form. The stronger radia- 
tive pumping populates higher rotational levels, which ac- 
counts for the bell shape of the curves. The same trend 
occurs at higher rotational levels, for which the curve max- 
ima are shifted towards higher densities. 

The observed excitation temperature of ~ 250 K is re- 
produced for the correct column densities, in a range of 
models, as shown on Fig. 1151 One typical example is illus- 
trated in Fig. [THl It can be clearly seen however that the 
J = 4 point lags under the J = 3 — 5 curve. This is directly 
linked to the imposed ortho-to-para ratio on the formation 
on grains. This is because the high radiation field results 
in a photodissociation time scale which is an order of mag- 
nitude larger than the conversion time scale from ortho-to- 
para H2 through reactive collisions, with either H or H + . 
The only way to reproduce the J = 4 column density is 
to assume a formation ortho-to-para ratio that is different 
from the statistical equilibrium value of 3. Fig. [TC] shows 
that an initial ratio of 1 provides robust results. We are 
however unaware of any theoretical or experimental result 
that would support such an initial ratio. 

We conclude that a pure steady-state PDR model is 
able to reproduce the excitation of H 2 at J > 2, observed 
towards HD 102065, with two ad hoc assumptions: the pres- 
ence of warm and dense gas close to the star at a pressure 
> 10 6 Kcm -3 , and an H2 ortho-to-para ratio, at formation, 
of 1. We checked the possibility that the warm H2 may be 
in a circumstellar disk. Upper limits on H2 column densities 
from the (3 Pictoris debris disk, are sever al orders of mag- 
nitud e lower than the HD 102065 values (jLecavelier et all 
2001). The warm H2 column densities observed towards 
HD 102065 are comparable to those detected from circum- 
stellar disks abo ut the young AeBe Her big stars HD 100546 
and HD 163296 (jLecavelier et alll2003l ). but these disks are 
traced by strong near to mid-infrared emission that is not 
ob served in HD 1 2065 

iGillmon et all (|2005f) and IWakkeH (|2006[ ) analyzed a 
large set of FUSE Galactic H2 detections obtained towards 
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Fig. 17. Fraction of H2 in the J=3 level versus the total H2 
column density. The data arc taken from observations of 
Galactic H2 towards extragalactic sources as analyzed by 
iGillmon et all (|2005l) . The three data points at N(H 2 ) > 
10 20 cm" 3 corre spond to the three Chamaeleon stars of 
iGrv et all (|2002l ) including HD 102065 (filled square). 

extragalactic sources. IWakkerl (|2006t ) showed that H 2 exci- 
tation, for column density ratios between the 4 first J levels, 
follows a systematic trend with increasing total H2 column 
density. HD 102 065, and the two remaining Chamaeleon 
stars studied by IGrv et alj (|2002D . fall on the same trend 
extending to yet higher N(H 2 ) values (Fig.[T7]). This agree- 
ment indicates that H2 excitation observed towards the 
Chamaeleon stars, fit with observations towards extragalac- 
tic sources and is thus unlikely to be due to matter heated 
by the stars. Qualitatively, the trend observed in Fig. [171 is 
in agreement with the fact that the gas temperature and 
the H2 radiative pumping increase, as the H2 column den- 
sity decreases. But, based on the HD 102065 results, we 
anticipate that a warm out-of-cquilibrium H2 component 
is required to quantitatively reproduce the data. We con- 
clude that the presence of some H2 at temperatures higher 
than the equilibrium temperature set by UV and cosmic- 
ray heating is a general characteristic of diffuse molecu- 
lar gas in the Solar Neighborhood. Modeling of the wide 
sample of FUSE Galactic H2 measurements, is required to 
statistically quantify the fraction of out-of-equilibrium H2 
gas. This warm molecular gas traces the local dissipation 
of turbulent kinetic energy. It is the dissipation of turbu- 
lence that creates the non-equilibrium chemistry discussed 
in Sect. 1531 

5. Conclusions 

We have gathered a wide set of observations characteriz- 
ing the diffuse molecular cloud, observed towards the star 
HD 102065. These observations provide independent con- 
straints on physical conditions, which are analyzed with the 
latest version of the Meudon PDR code. 

A single density (?ih = 80 cm" 3 ) slab of gas, bathed in a 
low radiation field, accounts for most observations (molec- 
ular fraction, gas temperature inferred from the N(H2, J = 
1)/N(H2, J = 0), cooling in the C + line, C 1 abundance and 
excitation). This model provides a physical and chemical 



reference upon which more elaborate interpretations can 
be developed to account for observables not reproduced by 
the model, namely the column densities of CH, CH + and 
H 2 in its excited (J > 2) levels. 

We consider the possibility that high J H2 is associated 
with matter close to the star. We place constraints on this 
possibility using IRAS data on dust emission. We find a so- 
lution where dense (nn > 10 4 cm -3 ) gas with a high pres- 
sure (p/k > 10 6 Kern" 3 ) would be located within 0.03 
pc of the star. This solution requires the presence of high 
pressure gas close to the star and that the H 2 ortho-to-para 
ratio at formation is 1. Such a departure from the statistical 
expectation is supported by no theoretical or experimental 
study. 

The H2 excitation observed towards HD 102065 fits with 
the general trend observed from FUSE Galactic H2 obser- 
vations towards extragalactic sources. We conclude that H 2 
excitation in the J > 2 levels observed towards HD 102065, 
is unlikely to be due to matter heated by the star, but it is 
more likely characteristic of H2 in the Solar Neighborhood 
diffuse ISM. Our work supports earlier studies that pro- 
posed that H 2 excitation in the J > 2 levels, traces the 
presence of warm H2, heated by the localized dissipation, 
in space and time, of turbulent kinetic energy within dif- 
fuse molecular clouds. The warm H 2 is physically associated 
with the bulk of the molecular gas, traced by the J=0 and 
1 H 2 absorptions, and is required to account for the CH + 
column density. In addition, the warm H2 could contribute 
to the CH abundance, and the inhomogeneity of the CO 
abundance, as indicated by the comparison of absorption 
and emission spectra. 

This paper outlines a framework for modeling the large 
number of Galactic H2 measurements derived from FUSE 
extragalactic observations, which cover 6 orders of magni- 
tudes in H2 column densities. In diffuse clouds, the time 
scale of H2 formation, a few 10 7 yr, is long compared to 
dynamical time scales. Observed values of the gas fraction 
in H2 are thus not expected to match steady-state values. 
Modeling of the available data will quantify how observed 
molecular fractions scatter about model values, and bet- 
ter constrain the H 2 formation rate than the present study. 
Modeling of these observations will, in addition, statisti- 
cally quantify the presence of warm H2 , heated by localized 
dissipation of kinetic energy. 
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